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1. Introduction 

PSR B1259—63/LS 2883 is a member of the small class of very high energy (VHE; E > 100 
GeV) /-ray binaries which comprises only five known objects. PSRB1259—63/LS 2883 is unique 
for being the only /-ray binary for which the compact object is clearly identified as a pulsar. This 
pulsar with a spin period of ~ 48 ms and a spin-down luminosity of ~ 8 x 10^^ erg s^' is moving in 
a highly eccentric {e = 0.87) orbit around a massive Be star with a period of Porb = 3.4 years (1237 
days) s@- The companion star LS 2883 has a luminosity of L* = 2.3 x 10^^ erg s ^ Because 
of the fast rotation of the star, it has an oblate shape with an equatorial radius of Peq = 9JRq and 
a polar radius of Ppoie = 8. IP©, which in turn leads to a strong gradient of the surface temperature 
from Peq ~ 27,500 K at the equator to Ppoie — 34,000 K at the poles [Q]. 

The Be star features an equatorial circumstellar disk, which is believed to be inclined with 
respect to the pulsar’s orbit (see e.g. ||3|]), so that the pulsar crosses the disk twice each orbit. The 
circumstellar disk of a Be star is a decretion disk with an enhanced stellar outflow formed around 
the star. As shown in optical interferometry observations, these disks are symmetrical with respect 
to the star’s rotation axis (see e.g. 0). Circumstellar disks generate excess infrared (IR) emission 
produced through free-free and free-bound radiation, providing an additional IR photon field to the 
blackbody flux from the optical star. 

The dense medium of the disk is believed to play an essential role in the resulting emission 
from the system. The position of the disk can be localized based on the disappearance of the pulsed 
radio emission from the pulsar. The observed radio emission far from periastron consists only of the 
pulsed component [Q, but closer to periastron, at about tp — 100 d (tp is the time of periastron), 
the intensity of the pulsed emission starts to decrease and completely disappears at about tp — 20 d. 
It then re-appears at around tp -|- 15 d. This eclipse of the pulsed emission is believed to be caused 
by the circumstellar disk. It is accompanied by an increase of the transient unpulsed radio flux 
beginning at ~ ?p — 30d and reaching its maximum at ~ fp — 10 d. This is followed by a decrease 
around the periastron passage and a second peak at about fp -|-20d A similar behavior is 

observed also for the unpulsed X-ray emission. Close to the periastron passage the X-ray emission 
features two peaks at around 20 days before and after periastron with flux levels 10-20 times 
higher than during apastron, and a decrease of the emission at the time of the periastron passage 
itself. The X-ray data is very similar from orbit to orbit repeating the shape of the light curve 
very well (0 and references therein). These peaks might be connected to the crossing of the disk 
environment. 

At TeV energies the source was observed by H.E.S.S. around four periastron passages in 2004 
[§], 2007 0, 2010 [|T^], and 2104 (results not available yet). The TeV emission from the source 
shows a variable behavior around the periastron passage. Although the exposure of the source is not 
sufficient to draw firm conclusions, the combined light curves from all three observing campaigns 
show a hint of two asymmetrical peaks before and after periastron, which roughly coincide in time 
with the peaks of the emission in the radio and X-ray bands, and a decrease of the flux at periastron. 
The shape of the light curve might be a result of the interaction of the pulsar with the circumstellar 
disk. 

Eirst observations at GeV energies conducted by Fermi-LKE took place around the 2010 pe¬ 
riastron passage revealing quite unexpected results [|rT|, [1^]. Fertni-LKE detected a low flux from 
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the source close to periastron with a subsequent dissapearance of the source after periastron fol¬ 
lowed by a sudden flare (with ~ 10 times the pre-periastron flux level) 30 days after periastron. 
This flare lasted for about 7 weeks without any obvious counterparts at other energy bands. The 
flare is shifted in time with respect to the post-periastron peak at other energy bands. The nature 
of the flare is still not understood, but several possible explanations have been discussed in the 
literature |[I^, [I^, [l^ |l 6 , EM- New observations around the 2014 periastron passage revealed 
the GeV flare at the similar orbital phase establishing the repetitive behaviour of this phenomenon 
[18]. However, new observations did not show any significant emission close to periastron [|T^. 

The circumstellar disk of the companion star plays a crucial role in the variability of the y-ray 
emission from the system. The dense disk photon field should significantly contribute to the target 
photons for inverse Compton scattering from the source, enhancing the observed TeV emission. 
The abrupt change of density at the pulsar’s entrance and escape from the disk should also cause 
a change of the shape of the shock between the pulsar wind and stellar environment, which might 
be the reason of such spectacular events as the GeV flare. The high density of seed photons should 
also increase the opacity for 77 -absorption, followed by electron-positron pair production, which 
in turn might scatter again the disk photon field, generating secondary 7 -rays. This cascade process 
may cause the re-emission of the TeV IC flux at lower GeV energies. Moreover, because of the 
deflection of electrons and positrons in the magnetic field, secondary 7 -rays can be re-emitted into 
randomized directions. Therefore, even if the primary VHE photon was emitted into the direction 
opposite to the line of sight to the observer, it can still contribute to the observable flux at lower GeV 
energies through the cascade emission. The pair cascade emission in binary systems caused by the 
interaction of the primary very high energy photons with the stellar photons was studied in detail 
for several cases (see e.g. [j^, and references therein). These studies showed that, although 
environments of the binary systems fulfill all the requirements for effective pair cascading, the 
resulting spectrum is in conflict with the GeV data observed by Fermi-LAI, since the expected 
cascade emission peaks in the Femii-LP^I energy band with very constraining upper limits. The 
77 -absorption of the TeV 7 -rays by stellar photons was also suggested as a possible explanation of 
the variability of the TeV flux across the orbit [^. The idea is that at periastron, when the pulsar is 
at the shortest distance to the star, the absorption should be the most effective, which would provide 
the decrease of the flux. However, showed that absorption only cannot explain the TeV light 


curve. 


2. Cascade Emission in the Circumstellar Disk 


We used the cascade Monte Carlo code of ||^, to calculate the angle-dependent cascade 
spectra for a variety of different input parameter sets within a parameter space motivated by the 
known properties of the Be star and its circumstellar disk in the PSR B1259—63/LS 2883 binary 
system. In these calculations we assume that the pulsar is a point-like source emitting isotropically 
in all directions. This assumption is reasonable if the wind termination shock at which leptons 
injected by the pulsar are efficiently accelerated and isotropized is close to the pulsar and the 
density of target seed photons is high enough to lead to IC scattering in the direct vicinity of 
the shock, which is the case for PSRB1259—63/LS 2883 [p4|]. We approximate the spectrum of 
the VHE 7 -rays as generated by IC scattering by electrons as a power law distribution with an 
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exponential cut-off. In our simulations we consider a mono-directional beam of photons to isolate 
all the geometrical effects. The source would emit the same beam of photons in every direction, 
and in the case of efficient cascading even those photons emitted in the opposite direction from the 
observer can contribute to the resulting observable spectrum. The orientation of the magnetic field 
would be differenf for primary photons emitted in differenf directions. The disk is approximated 


by a grey-body wifh the energy density Ud and temperature Td = 0.6 [25, M, 27] and the shape 


of the disk is approximated by a cuboid whose side lengths are chosen in a way to correspond to 
the real size of the disk. All the details on model assumptions and their justification can be found 
in 0. 

The dependence of the cascade emission on the viewing angle, magnetic field sfrengfh and 
orienfafion, energy densify in fhe disk, and source location within the disk were discussed in detail 
in [Mj. It was also shown that the GeV flare is very unlikely fo be genereafed by fhe cascade emis¬ 
sion. However, cascade emission mighf sfill confribute fo fhe fainf GeV emission observed close 
fo periasfron during fhe periasfron passage in 2010. The specfrum presenfed by [l^j is consisfenf 
with a substantial contribution from cascade emission. The energy of the peak in the spectrum of 
the cascade /-rays emitted in the forward direction only slightly depends on the energy density and 
magnetic field [M] and is in fhe energy range of 15 — 30 GeV. Figure |T] shows fhe comparison of 
fhe GeV-TeV dafa wifh fhe simulafion resulfs of fhe cascade Compton radiation from fhe primary 
phofon beam direcfed towards fhe observer emitted in fhe forward direcfion info a cone wifh an 
opening angle 11° (0.98 < ft < 1). This corresponds fo fhe cascade emission radiated in fhe di¬ 
recfion fo fhe observer from fhe primary /-rays emitted wifhin fhe cone of fhe same opening angle 
assuming fhaf fhe magnefic field affecls fhe cascades initiated by every primary photon idenfically. 
This assumption is valid for sufficienlly small values of fhe opening angle, such fhaf fhe orienfafion 
of fhe magnefic field wifh respecf fo the direction of the primary photon does not change much and 
its influence on fhe shape of fhe cascade emission specfrum is negligible. The simulafion was per¬ 
formed for Ud = 20 erg/cm^, magnefic field B = 0.1 G wifh j3 = B^/By = 10 (B lies in fhe xy-plane). 
We decided fo chose a magnefic-field direcfion closely aligned wifh fhe primary photon beam (and 
fhus, fhe direcfion fo the observer). In this way, the magnetic field has very little influence on fhe 
cascade developmenf and on the resulting Compton emission. Otherwise, the cascade radiation 
emitted in other directions can become higher than the forward emission and a cone with a bigger 
opening angle should be considered, which in turn requires an accurate calculation of the magnetic 
field orienfafion for every primary phofon. 

However, fhe specfrum of fhe GeV emission during fhe "low-flux" period derived by an inde- 
pendenf group feafures a completely differenf specfral shape [ITTl]. If reveals a significanl flux af 
0.1 — 1 GeV wifh no emission above 1 GeV. In fhis case fhe specfrum rejecfs fhe possibilify of a 
significanf cascade confribufion, because even a very small level of absorpfion of fhe TeV /-rays 
and re-emission af lower energies would violafe fhe Fenni-LKP upper limif in 1 — 100 GeV energy 
band. The Fermi-LKF upper limif in fhis case consfrains fhe energy densify in fhe disk fo uj < 8 
erg/cm^ [Mj. Unfortunately, no significanf emission close to periasfron was defected during fhe 
2014 periasfron passage [p^], fhus, fhe specfral discrepancy remains unresolved. 

If should be nofed fhough fhaf fhe simulations presenfed above are performed considering only 
fhe disk radiation field and neglecting fhe stellar radiafion field. Therefore, fhe comparison of fhe 
model wifh fhe "periasfron" emission presenfed on Fig. |T] should be freated only as a hinf fhaf fhe 
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Figure 1: Comparison with the data. The black dashed line represents the intrinsic 7 -ray spectrum. The 
blue solid line represents the spectrum of the forward emission (0.98 < /i < 1) modihed by gamma-gamma 
absorption and cascading. Assumed parameters: B = 0.1 G, j3 = 10, Ud = 20 erg/cm^. Black squares repre¬ 
sent the 2011 H.E.S.S. data [ |l^ and red triangles show the GeV data close to the 2010 periastron passage 
as reported by [p^ . 


GeV emission may be explained by the cascade radiation (even if we consider only disk photons). 
This implies that taking into account stellar photons might probably provide a better fit of the data 
and decrease the estimate (upper limit) of the disk energy density. A detailed study of the cascade 
development in the system which would properly account for all relevant effects is in progress. 


3. Impact of the 77 -Absorption on the Observed TeV Light Curve 


The absorption of the 7 -ray emission caused by the interaction with the disk radiation field 
might have a significant impact at the observed TeV light curve. The resulting integrated flux 


above 380GeV after absorption caused by the stellar radiation field was calculated in []21p. It 
was shown that the absorption due to the stellar radiation plays only a minor role in the observed 
variability. However, the disk radiation was not considered in those calculations. Using the estimate 


of the resulting integrated flux after the absorption on stellar photons only, F^, calculated in [ 21 ] we 
calculated the resulting flux after the total absorption due to both stellar and disk radiation fields as 


F{E > 380GeV) ~ F,{E > 380GeV)e~^’T' = E,{E > 380GeV). 


D ^7, mm / 


^7, max ^7, min 


(3.1) 
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assuming the optical depth of the disk radiation field d(£) only weekly depends on energy 
(£y = Ey /is the normalised /-ray energy) in the energy range Sy.min = 380 GeV and Ey^^^x = 
lOTeV and can be replaced by a constant averaged value Zyy. For these calculations we used the 
same orbital parameters as used in [ ^ ] : orbital period ■Porb = 1236.7, eccentricity e = 0.87, perias- 
tron longitude ft) = 138.7°, and inclination angle / = 35°. The disk inclination angle is assumed to 
be id = 10°. The disk is assumed to be perpendicular to the major axis of the orbit with a constant 
width of 10*^ cm (corresponds to 1° half-opening angle at a distance of 457?*) and constant energy 
density of 8 erg cm^^ (the highest energy density for which the Eermi-LAJi upper limits are not 
violated, see previous subsection). The expected light curve after the total absorption assuming 
a constant initial /-ray flux is shown in Fig. ^ by the dashed line^ The solid line represents the 
flux after absorption taking into account only stellar photons as reported by pT|]. However, the 
assumption of the constant optical depth underestimates the effect of the //-absorption in the disk 
as the optical depth d(£) peaks at ~ 300 GeV. The resulting light curve after absorption in the 
disk only taking into account dependence of the optical depth on energy is shown by the dotted line 
in Fig. II 

It can be seen that the circumstellar disk might play a crucial role in the flux variability around 
the periastron passage and that the observed flux variability can be explained by the absorption. 
The time shift of the observed preperiastron peak with respect to the modelled light curve might be 
explained by the rotation of the disk normal with respect to the major axis. The assumption of the 
constant width of the disk results in an overestimation of the absorption close to periastron, since 
the disk should be thinner close to the star. However, the energy density in the disk is expected to 
be highest close to the star. The combined effect might result in a shallower minimum near perias¬ 
tron than predicted by our calculations. The light curve has been modelled under the assumption 
that the primary /-ray flux is constant. This is a valid approximation only in the case when IC 
scattering occurs in the saturation regime which is believed to be true around periastron. Farther 
from periastron the primary flux is expected to decrease. In the case when IC scattering is not in 
the saturation regime, local maxima are expected at the disk crossings and at the periastron. This 
would completely change the expected light curve. With the current quality of the observed TeV 
light curve it is difficult to draw any final conclusions. New observations conducted in 2014 might 
bring some clarity and give answers to some of these open questions. 

A critical test of the absorption model could be derived from the study of the spectral variability 
of the TeV emission. The TeV spectrum is expected to soften when the pulsar is moving through 
the disk for the first time, then the spectral shape should remain almost unchanged with a slight 
variability due to the absorption in the stellar radiation field until the pulsar moves out of the 
disk after periastron. After that the spectrum should harden again. Unfortunately, the VHE /-ray 
observations during previous periastron passages did not provide enough statistics to draw firm 
conclusions about the spectral variability. Although there are some indications of spectral changes 
both during the 2004 and the 2007 |H| periastron passages, these changes are not significant. 
However, if the spectral variability is real on both occasions, it suggests a hardening of the spectrum 
towards periastron, which is in contradiction with the effect we expect in the case of effective 

'Note, that there was a mistake in calculation of the angle between the disk normal and line of sight in the original 
paper p^]. Corrected results presented here show a slightly stronger effect of the disk absoprtion. 
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Figure 2: VHE 7 -ray flux modulation caused by 77 -absorption in PSRB1259—63/LS 2883. The solid line 
represents the normalized integrated flux above 380 GeV after absorption caused by the stellar radiation field 
as calculated in The dashed line shows the same flux taking into account also the absorption caused by 
the disk radiation field assuming the 77 optical depth is constant above 380 GeV. The dotted line represents 
the normalized integrated flux above 380 GeV after absorption caused by the circumstellar disk only taking 
into account the dependancy of the 77 optical depth on energy. See the main text for details on the geometry 
and other parameters. The data points represent the combined H.E.S.S. light curve from three observed 
periastron passages as reported in [|^. The integrated flux above 1 TeV is extrapolated down to 380 GeV 
using a photon index of 2.8. The same flux normalisation factor as in [^] of 10^*^ cm^^ s^* is used. 


cascading. 

The study of the 77 -absorption in the disk which would take into account a proper geometry 
of the disk and a realistic distribution of IR photons within the disk corresponding to the observed 
IR emission from the star is currently in progress. 
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